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A NOTE ON THE OPTIMUM AVERAGING TIME OF WIND 
INFORMATION FOR CONVENTIONAL AIRCRAFT LANDINGS 


By W. R. SPARKS and BARBARA KEDDIE 


. From detailed analysis of wind measured by a single anemometer 10 metres 
above open level ground, and consideration of Shiotani’s observations with a horizontal array 
of anemometers, it is concluded that the optimum averaging time for wind measured by a 
single anemometer should not be less than four minutes. 


Introduction. At some stage during an aircraft approach it is necessary 
to pass to the pilot the wind expected in the touch-down area at the time 
his aircraft will reach that area. Since it is not possible to site a conventional 
anemometer in the touch-down area this necessitates the extrapolation of 
measured wind information in both time and space, and prompts the question : 
what can be measured with the normal airfield anemometer to give the best 
estimate of the wind that an aircraft will encounter in a different position 
and at a later time? 


Before it is possible to begin to answer that question, what information 
the pilot requires must be known. There is no general agreement on this, 
but discussions with aviation experts suggest that the average wind over a 
period of 4 or 5 seconds just before touch-down when the aircraft is descending 
from about 30 m to 15 m is of critical importance. During that period the 
aircraft would travel about 300 m. If the validity of Taylor’s hypothesis of 
frozen turbulence is assumed (Panofsky et alii), the analogous measurement 
from a stationary anemometer is the average speed of a 300-m run of wind 
(Uso0)- 

Let the wind measured at the anemometer be averaged over a time T 
to obtain a value U7 and that average be used as an estimate of Ugo. at some 
later time. How wrong can it be? A good measure of the errors is the root- 
mean-square difference between Ugo) and Ur divided by ogo9, the standard 
deviation of Uso. over VV observations, i.e. 


(1/¢390)V= [(Uso0—U7r)*/N]. 


Results and discussion. From observations made with an anemometer 
10 m above flat open country at Cardington (near Bedford), the above ratio 
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has been calculated for averaging times 7 from 4 minute to 15 minutes. 
(The observations are described in Appendix I.) The ratio does of course 
depend on the interval between the observation of Ur and the time that 
Uso9 is required, see Appendix II. Two intervals have been considered, 
10 minutes and zero minutes. The results of these calculations are shown 
n Figure 1. 





Averaging time, T \minutes 


FIGURE I—ERRORS IN ESTIMATING Uso9 FROM U7 (AVERAGED OVER SIX DATA 
RUNS) 


A For a 10-minute interval between Uy and Ugo 

B_ For a zero time interval between U7 and U. 

C_ Theoretical curve for zero correlation between short-period means 
at two points 


Curve A, which gives the results when there is a 10-minute interval between 
the observation and its use, i.e. between the end of the averaging time 7 
and the beginning of the run of wind, shows that the minimum error is made 
when an averaging time of 4 to 5 minutes is used. Curve B shows that even 
when the information is used immediately after the observation, no increased 
accuracy is obtained by reducing the averaging period to any less than the 
4 minutes suggested by curve A. 


Curves A and B show the results that would be obtained if the anemometer 
could be sited in the touch-down area. Consider now the effect of measuring 
some distance from that area. The reduction in accuracy is a function of 
the correlation between the wind at the measuring point and the wind in 
the area of interest. This correlation depends on the separation between 
the instrument site and the area of interest and the averaging periods used. 


On a normal airfield the anemometer would be several hundred metres 
from the touch-down area. The results from Shiotani’s* observations with 
a horizontal row of anemometers suggest that the correlation between the 
airfield anemometer’s indication of short-period means and simultaneous 
observations in the touch-down area could well be zero. 


Curve C, determined as in Appendix II, shows the errors in estimating 
Uso9 in the touch-down area if the long-period means (say hourly running 
means) at the anemometer and in the touch-down area are equal and constant 
but the variations of the short-period means about the long-period mean 





Meteorological Magazine, 100, 1971 131 


are uncorrelated at the two points. The variances of Us99 and Uy used to 
construct curve C were the values measured at Cardington during Run No.3 
(see Appendix I). 


If the anemometer is sited a considerable distance from the touch-down 
area, curve C is applicable and it can be seen that any attempt to use an 
averaging period of less than the 4 minutes suggested by curve A would 
increase the errors even if the information were used immediately after the 
observation. 


It has been demonstrated that, even if the wind information required 
by the pilot is a mean over a shorter period, the best estimate of the required 
information that can be obtained from a conventional anemometer is a mean 
value over a period of at least 4 minutes. It would also be possible, and may 
be desirable, to give the pilot some measure of the uncertainty of the estimate. 


Although these conclusions are based on limited data (Appendix I) they 
are in excellent agreement with those of Rijkoort and Wieringa® and it is 
reasonable to assume that they are applicable to the well-stirred conditions 
in neutral equilibrium which accompany strong winds near the ground. 
The sudden arrival of squalls associated with the near approach of thunder- 
storms or the passage of a cold front would probably prove to be a source 
of large errors, but such errors would be expected in any method relying on 
observations from a single anemometer. There is however, a need for 
information on the effects of spatial separation between the anemometer site 
and the touch-down area. 


Appendix I 


Description of data analysed. The data were obtained from a three-cup 
anemometer with a distance constant of 5 m and a wind vane with a sine/ 
cosine potentiometer. The data analysed were the westerly components of 
winds from directions between 240 degrees and 280 degrees. The anemometer 
was sited 10 m above open level country at Cardington. The site was fully 
described by Giblett* in 1932 and has not changed significantly since that 
time. The data were logged on 31 March 1969. During most of the period 
Cardington was in a warm sector and the winds were sufficiently strong to 
ensure near neutral stability. The mean wind was not, however, steady; 
during one data run it changed by about ro kt in three minutes. 


The data were logged at one-second intervals with each data run lasting 
about 50 minutes. The main features of the wind during each run are shown 
in Figure 2 by plots of consecutive one-minute means for Runs 2, 3 and 5. 
The features of Runs 1, 4 and 6 are similar to those of Run 3. 


Figure 1 (page 130) shows average values over all six runs of the ratio 
(1/¢390)V =[(Uso0.—Ur)?/N]. 


Figures 3 and 4 show the same information for individual runs. 
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FIGURE 4—ERRORS IN ESTIMATING Uso) FROM Uy FOR A ZERO TIME INTERVAL 
BETWEEN THE TWO OBSERVATIONS 


Appendix II 
The construction of the theoretical curve for zero correlation 
between short-period means at two points. Uy, can be expressed as 
the sum of a long-term mean Uy appropriate to the mid-point of the interval 
T and the departure U’7 from that mean, i.e. Up = Uy + U7. 


Also Uso. can be expressed as the sum of a long-term mean Uj, appropriate 
to the mid-point of the interval of time for Uy, and the departure Uo, 
from that mean, i.e. Us99 = Um + U’ soo: 


Consider the term 
I T 
F= . s 2 (Us09 — Ur)?, 
where WN is the number of observations of U399; 


I 
“ N 2 (Um — Um + U'g99 — U'r)? 


I ay ts 
= —— 2( Um — Um)? + => (Um — Um) (U'g99 — ;U'r) + 
N N 


I 
4. wv BU s00 — U'r)*. 


The second term will be zero if the difference between the two departure 
terms is not correlated with the trend in the long-term mean. 


I I 
Then F = —-E(Um — Um)* + Se E(U'soo)* + Ge E(U'n)* — 


I 
om “N_ x2 U’' s00 U'r 


I 
Be 2 X (Um—Um)* + 6* 399 +0? 7 —21rGg9997, 
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where r is the correlation between the variations of the short-period means 
from the long-period means. 


This final equation for F may be used to explain the construction of curve 
C in Figure 1. If the long-period means at the anemometer and in the touch- 
down area are equal and constant, i.e. Um=Uy=constant, but the variations 
of the short-period means about the long-period means are uncorrelated at 
the two points, i.e. r=o, the expression for F becomes 


I 
F = 6%399 + O87 = WV 2 (Us00— Ur)? 


Hence in the construction of curve C, the ratio 


(1/300) V 2 [(Us00 — Ur)*/N] 


becomes 


V (1 + o7*/og90*)- 
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SOME ASPECTS OF WIND INFORMATION REQUIRED IN THE 
LANDING OF AIRCRAFT 


By BARBARA KEDDIE 


Summary. Wind measurements were made each second with a single anemometer over a 
period of about an hour for six independent wind régimes. Detailed analysis was made of 
the maximum difference between a mean speed Up measured over 44 minutes and the speed 
Ugo9 10 minutes later over a 300-m run of wind. The maximum value of Ug9)— U7 in various 
samples of sizes from 4 to 20 minutes was predicted by assuming a normal distribution. 
Theoretical curves were derived which were in reasonable agreement with curves drawn from 
the analysed data. Actual values were also compared with predicted values of (i) maximum 
deviation of Usg)— Ur from the mean in various samples, and of (ii) the maximum difference 
in various samples. 

Finally comparison was made of data from different sampling periods or runs of wind so 
that if the maximum departure of one sample from the mean wind was known the maximum 
departure for other samples could be estimated. 


Introduction. Pilots and others concerned with the operation of aircraft 
need a description of the airfield low-level wind which is used as a forecast 
of the wind expected a few minutes later as an aircraft touches down. At 
present the pilot is given the ‘surface wind’ by the controller about 10 minutes 
before touch-down. The problem is to find what information could be given 
to the pilot which would best describe conditions 10 minutes later. 
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The behaviour of the wind in the last 30 m of descent, in particular between 
30 m and 15 m, is most important to a descending aircraft whether a manual 
or an automatic landing is being carried out. An aircraft normally travels 
about 300 m while descending from 30 m to 15 m, thus the 300-m run of 
wind defines one scale of eddy of importance to aviation. Therefore it 
would be useful to find a measurement which is representative of the 300-m 
run of wind 10 minutes later. 


It has been assumed here that Taylor’s hypothesis of frozen turbulence 
holds. According to this hypothesis a space correlation in the direction of 
mean wind, x, can be determined from the time correlation function by 
replacing time t by x/u, where u is the mean wind speed (Panofsky'). 


By comparing the mean speed over a 300-m run of wind with that measured 
over a period T ending 10 minutes earlier, for values of T between 14 and 
15 minutes, the optimum averaging time T has been found to be at least 
4 to 5 minutes (Sparks and Keddie?). For the limited data considered here, 
the best estimate of the mean wind speed over 300 m to be expected in 10 
minutes’ time was taken to be the mean speed over the last 44 minutes, as 
longer averaging times did not appear to give an improvement in the accuracy 
of the estimate. 


A pilot also requires information on the gustiness of the wind, but the 
anemometers in general use on airfields cannot measure gusts on a very small 
linear time-scale. At Cardington there is an anemometer placed 10 m above 
flat open country from which a read-out every second can be obtained (see 
Sparks and Keddie, Appendix I). By adding together the run of wind for 
successive seconds until a 300-m run has been obtained and dividing by the 
number of seconds, the mean wind speed over a 300-m run of wind can be 
found. Similarly mean speeds can be found for other lengths of run of wind. 
By comparing the wind extremes for different runs of wind a relationship 
might be determined whereby, if the maximum speed of one run of wind 
was measured, that of a different run of wind could be estimated. 


Data available. The data used were obtained during six Runs at 
Cardington on 31 March 1969; details are given in Table I. A warm front 
passed through the station just after 08 cmt and the associated cold front 
passed between 17 and 18 omrt followed by a trough soon after 19 GMT, 
i.e. near the start of Run 5. The data consist of a read-out every second of 
the westerly component for winds from between 240 and 280 degrees, 
measured in metres per second at a height of 10 m (see Sparks and Keddie, 
Appendix I). 


TABLE I—DATA OBTAINED DURING SIX RUNS AT CARDINGTON ON 3I MARCH 1969 


Run number Start time Duration Mean speed Mean direction 
GMT min m|s degrees true 
1313 48 13°7 240 
1555 55 12°0 240 
52 11-2 245 
52 10°5 245 
53 11-7 ato" 
52 11-2 265 
* Discontinuous, 250° at first becoming 270° after 3 minutes gradually backing to 255° by 
1945 GMT. 





136 Meteorological Magazine, 100, 1971 


For successive 300-metre runs of wind the mean speed Ugg9 was found 
and compared with U7, the mean wind speed over a 44-minute period ending 
10 minutes earlier. The mean of the differences (U399—U7r)m, and standard 
deviation of the differences og99-r, were calculated for each Run. If the 
differences are normally distributed a prediction within given confidence 
limits can be made of the maximum (U39,—U7r) likely in a given period 
with a steady régime. 


Comparison of values predicted for the maximum (U,,,.— Ur), 
assuming a normal distribution, with those actually obtained. 
Assuming a normal distribution, the probable maximum value of (U39,.— Ur) 
in a sample of V values may be defined either as (i) the value likely to be 
equalled or exceeded once in a sample of WV (Durst*), or as (ii) the value 
likely to be exceeded once in a sample of 2N (Brooks and Carruthers‘). 


Samples of 8, 16, 24 and 32 values of (Us9,—Uz7) were considered, 30 
random samples of each size being taken for each Run. The mean maximum 
of (Us9.—Ur) was found for each set of 30 samples and standardized by 
subtracting (Us99—Uz7)m and dividing by o399-7 to facilitate comparisons 
with predictions based on normal distributions. Figure 1 shows the results 
together with curves, normal (1) and normal (2), determined from the normal 
distribution for definitions (i) and (ii) of the probable maximum value. 
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FIGURE I1—COMPARISON OF MEAN MAXIMUM (U39,—U7) OBTAINED IN RANDOM 


SAMPLES OF VALUES OF (U39,.—U7) wiTH THE MAXIMUM PREDICTED ASSUMING 
A NORMAL DISTRIBUTION 





Definition (i) is perhaps more appropriate to the way in which the maximum 
is found from the samples. The curves are seen to be in reasonable agreement, 
except for Runs 2 and 5. Near the beginning of Run 5, there was a large 
increase in wind speed, about 10 kt in 3 minutes, and therefore there were 
some values of (U39,—U7) at the beginning of the Run which were much 
larger than the mean. The probability of one of these large values being 
included in a random sample is quite high, especially in the larger samples, 





Meteorological Magazine, 100, 1971 137 


and so the ‘mean maximum’ is likely to be high. On the other hand, in Run 2 
there was a decrease in wind speed in the middle of the Run, leading to 
several large negative values of (U39,—U7). This lowered the value of the 
mean maximum and also led to a large standard deviation, and hence the 
values plotted for Run 2 seem abnormally low. 

Thus it seems reasonable to assume a normal distribution, at least in cases 
where there are no rapid and permanent changes of the wind, and to use 
definition (i). However, in practice it is not random values of (U39.—Ur) 
that one is concerned with, but the maximum gust which is likely in a certain 
length of time, or in a certain number of consecutive values of (U399—U 7). 
Since a 300-m run of wind took from 30 to 40 seconds, sample sizes of 8 to 
32 consecutive runs of wind represent time intervals of from about 4 to 20 
minutes, which is adequate for the consideration of prediction for 10 minutes 
ahead. Figure 2 shows the results obtained using samples of 8, 16, 24 and 32 
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FIGURE 2—COMPARISON OF MEAN MAXIMUM (U39,—U7) OBTAINED IN SAMPLES 
OF CONSECUTIVE VALUES OF (U3 9,—U7) wiTH THE MAXIMUM PREDICTED 
ASSUMING A NORMAL DISTRIBUTION 


consecutive values of (U39,—U7). The number of samples of each size varied 
as the Runs were not all of the same length; for the shortest Run there were 
23 samples of 8, 19 of 16, 15 of 24, and 11 of 32; and for the longest Run there 
were 34 samples of 8, 30 of 16, 26 of 24, and 22 of 32. The curves are not in 
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very close agreement with either of the normal curves, although they show a 
similar shape. The lack of agreement could be because of changes in the 
mean wind speed over the Run. To overcome this effect the range was 
considered instead of just the maximum value, i.e. for each sample the 
difference between the maximum and minimum values of (Us99—Ur) was 
found and the mean taken over a number of samples. By dividing this mean 
range by twice the standard deviation, ¢s99-7, standardized curves can be 
obtained (Figure 3) which can be compared with the normal curves used in 
Figures 1 and 2. 


2:2 


~—Normal(? 


Normal(1) 


Run3 








I 1. 

16 24 

Number in sample 
FIGURE 3—COMPARISON OF MEAN RANGE OBTAINED IN SAMPLES OF CONSECUTIVE 
VALUES OF (Ug99— U7) WITH THAT PREDICTED ASSUMING A NORMAL DISTRIBUTION 


The standardized curves obtained in Runs 1, 3 and 4 are quite close to 
the normal (1), and these were Runs where conditions were fairly steady. 
Run 5 is again abnormal, the very low values being due to the large standard 
deviation. 


Although the values in Figures 2 and 3 are in some cases rather far from 
the normal, the slopes of the curves appear to be similar, implying that the 
ratios of points on the curves are similar for each Run. Figures 4 and 5 show 
the ratios of values for samples of 16, 24, and 32 to those for samples of 8, 
corresponding to Figures 2 and 3 respectively, for each Run and for the 
normal curve (1). The agreement between all the Runs and the normal in 
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FIGURE 4—RATIO OF MAXIMUM DEPARTURE OF (U3 9,—U7) FROM THE MEAN 
FOR DIFFERENT SIZED SAMPLES TO THAT FOR A SAMPLE OF 8 
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FIGURE 5—RATIO OF MEAN RANGE FOR SAMPLES OF (U3)9,—U7) TO THAT FOR 
A SAMPLE OF 8 


Figure 4 is very good. The agreement between the normal and Runs 1, 3 
and 4, when conditions were fairly steady, in Figure 5 (which considers range 
instead of maximum) is extremely good, but for the remaining Runs the result is 
not as good as in Figure 4, indicating that for general use the range is not such 
a good parameter as the maximum value as far as ratios are concerned. This 
shows that if the maximum gust in a sample of a certain size is measured then 
the maximum likely in larger samples can be estimated by assuming a normal 
distribution and using the appropriate ratio. Almost equivalent to this is 
the statement that if the maximum gust in a certain period of time is 
measured then the maximum likely in a longer period can be estimated. 


Table II(a) gives an indication of the errors obtained in using the normal 
assumption in practice. For each Run, the maximum positive deviation of 
(Us99—Ur) from the mean was found for consecutive samples of 8 and for 
the larger samples, each of these starting with the sample of 8. The values 
obtained were compared with those predicted — using the ratios shown in the 
normal curve in Figure 4 — from the data of the sample of 8. The table 
shows that, apart from some large errors in Run 5, the errors in all other 
forecasts are within the range -+-2 m/s. 
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TABLE 1(a)—ERRORS IN THE MAXIMUM DEVIATION OF Ugo99—Ur FROM THE 
MEAN IN VARIOUS SAMPLES AS PREDICTED* FROM THE DATA OF A SAMPLE OF 8 
Maximum deviation of Ugo9— Ur Predicted maximum deviation 
Run from the mean minus actual 
number Sample of Sample of Sample of Sample of Sample of Sample of Sample of 
8 16 24 32 16 24 32 
metres per second metres per second 
1°17 1°17 It] 117 0°39 0°59 0-73 
—0o'60 1:07 1°07 1*50 —0°27 —O'17 —0°53 
1°07 1°07 1*50 1*50 0°35 Or! 0°23 
0°62 1*50 1*50 1°50 —o-68 —0°57 —0"50 
210 2°10 2°10 2°10 0-69 1°05 1°30 
2:03 2°03 2°03 2°03 0°67 102 1°26 
08. 08. 084 08. 0:28 0°42 0°52 
o'5 05 o-7 1°3 0-18 0-10 —0°45 
—0'25 0°74 1° 1-69 —0-4! —o-98 — 1-28 
0°20 0°46 2°03 2°03 —O-1g —1'73 —1°7! 
0-46 2°03 2°03 2°03 —1.42 — 1°34 —1°28 
2°03 2°03 2°03 2°03 0°67 1°02 
0:98 0-98 0-98 Ivll 0°32 0°49 
O15 I'll 113 1°49 —o'9g! —0-g0 
Ind 1°13 1°49 1°49 0°35 o18 
1°13 1°49 1°49 1°49 ool O21 
4°83 4°83 4°83 4°83 1°59 2°42 
4°47 4°47 4°47 4°47 1°48 2°24 
1°31 1°31 1°31 1°31 0°43 0-66 
0-4! o-4! o-4! o'4!1 O14 O°21 
—0°53 0°34 0°34 0°34 0-36 0°46 
0°34 0°34 0°34 0°44 Orr O17 
2°79 2°79 2°79 2°79 0°92 1°40 
0°97 0°97 1°58 1°58 0°32 —O'12 
0-28 1°58 1°58 1°58 —121 —1°16 
1+58 1°58 1°58 1°58 0°52 0-79 
0°83 1°33 1°33 1°33 —0°23 —0-08 
* By multiplying by a factor based on the probable maximum values in a sample size if 
distribution is normal. 


(Uso90.—Ur)m is the mean value of (U3,,— U7) over the whole Run which 
it is to be hoped will be approximately zero. Table I1(b) shows the errors 
that would occur if the mean were assumed to be zero, i.e. just considering 
the maximum value in a sample and not subtracting (Us9,—Uzr)m. The 


table shows that this compares favourably on the whole with the first method 
on these occasions. 


Comparison of different sampling periods (distances). The 
Cardington data consist of a read-out every second, each being in fact a 
mean over about 0-9 seconds. The first 45 minutes of each Run were divided 
into a number of consecutive periods of T seconds duration, and the mean 
wind V7 in each T-second period was found. The maximum departures 
from V7 of the 1-second values and means of runs of wind over 40, 100, 200, 
300, and 400 m were found for each period T. By dividing each run of wind 
by the mean wind over the whole 45 minutes, the approximate equivalent 
number of seconds, t, was obtained for each run-of-wind distance. 


The mean wind speed over ¢ seconds may be written 


Vie =Vr+, 


where V7 is, as defined above, the mean wind speed over a longer period T, 
and vy; is the mean deviation from V7 over the period ¢ seconds. Let vg (max) 





Meteorological Magazine, 100, 1971 141 


TABLE 11(6)—ERRORS IN THE MAXIMUM DIFFERENCE Us99—Uzy IN VARIOUS 
SAMPLES AS PREDICTED* FROM THE DATA FOR A SAMPLE OF 8 


Maximum 
Run difference Predicted maximum difference minus actual 
number Sample of 8 Sample of 16 Sample of 24 Sample of 32 
metres per second 
1 1°43 O47 0-72 0°89 
—0°34 —o- —o-82 1°21 
13 0°44 0°24 0°39 
0" —0°59 — 0°44 —OSS 
1°74 0°57 0°87 1-08 
167 0°55 0°84 1°03 
0°48 016 0°2 0°30 
0°20 0°07 ~oa _ Fe 
—o'61 0°43 —0o-08 —0°34 
0°05 —0°24 — 1-80 —1°80 
7 —1°47 —141 —1-38 
1°88 0-62 0°94 1°17 
0°83 0°27 0°42 0°37 
—0'25 —o'38 —0°35 a 
o'71 or21 —0'02 o 
0°73 —o'12 or10 0°09 
5°40 1-78 2°70 3°35 
5°04 1°66 2°52 3°12 
1-88 0-62 0°94 1°17 
098 0°32 0°49 0-61 
0°04 —o'86 —o°85 —o8 
og! 0°30 0°46 O-4 
2°59 085 1-30 1°61 
0°77 0°25 —0'22 —0'13 
0°08 —1'27 — 1-26 —1'2 
1°38 0°46 o- o 
0°63 —0'29 a —Oll 
* By multiplying by a factor based on the probable maximum values in a sample size if 
distribution is normal. 


be the maximum value of % during a sample duration of T seconds. To 
compare different sizes of gusts, the ratio R(T, t: 6, s) may be used. This 
is defined as 


v7,t (max) 
Ug,s (max) ; 


R(T, t: 8,5) = 


where 6, s are respectively sample durations and averaging times, which may 
not be the same as 7 and ¢ (Brook and Spillane‘). 


In the Runs considered 100 metres is approximately equivalent to 12} 
seconds for the wind speeds, so the ratio R(T, t: T, 124) was calculated for 
each period 7, where ¢ is the number of seconds equivalent to a run of wind. 
As the 1-second value is actually a mean over about o-g seconds, the value of 
t in this case was taken to be 0-9. Two values of T were considered, 44 minutes 
and g minutes, and for each Run there were ten 44-minute periods and five 
g-minute periods. The mean values of the ratios for J = 270 seconds and 
T = 540 seconds were found for each value of t and for each Run, and the 
results compared with those obtained by Brook and Spillane which were 
worked out on the basis of ‘t-second gusts’ and not runs of wind (Figures 6 
and 7). The Brook and Spillane points in Figures 6 and 7 were obtained 
from Figure 4c in their paper, which gives values of the ratio R(T, t: T, 5). 
Two graphs were drawn, one for R(270, ¢: 270, 5) and one for R(540, t: 540, 5), 
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FIGURE 6—RATIOS OF MAXIMUM DEPARTURES FROM THE 44-MINUTE MEAN WIND 
OF MEANS OVER f/f SECONDS TO THE MEAN OVER A I00-METRE RUN OF WIND 
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and the two values of R( 7, 124: T, 5) obtained from the graphs. By dividing 
the Brook and Spillane values of R(T, t: T, 5) by the obtained value of 
R(T, 12: T, 5) the required ratios R(T, t: T, 124) were obtained, and these 
are plotted in Figures 6 and 7. The agreement in both cases is very good, 
the g-minute results being perhaps slightly better than the others. However, 
in both cases, the ratios obtained from the Cardington data are greater than 
those of Brook and Spillane for values of ¢ of 1 second and 5 seconds. This is 
almost certainly due to the presence of high-frequency ‘noise’ in the Cardington 
data. Thus, if the maximum departure of one sample from the mean wind 
over, say, 9 minutes can be measured, the maximum departure of other 
sample lengths can be estimated. 


Conclusion. A pilot can be given information which will give him an 
estimate of the mean speed and variability of the wind 10 minutes later. A 
mean wind speed measured over 44 minutes gives a reasonable estimate of 
the mean speed over a 300-m run of wind 10 minutes later. If the maximum 
departure of (U3 9,— U7) from the mean were measured over a short period, 
the maximum likely in a longer period could be estimated, assuming a normal 
distribution and using the appropriate ratio; for example, if the maximum in 
a sample of 8, i.e. about 4 minutes, were measured, then multiplication by a 
factor 1-62 would give an estimate of the maximum likely to occur in a sample 
of 32, i.e. about 16 minutes which would include the time for the sample of 
8 and the 10-minute period before the aircraft reaches the critical position, 
as well as the actual period of descent from 30 m to 15 m and then on to 
touch-down. The maximum gust appropriate to different aircraft sampling 
could also be estimated, using the methods of Brook and Spillane. 


The results outlined above are based on one day’s observations only. 
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EXTREME WIND SPEEDS IN THE COMMONWEALTH CARIBBEAN 
By H. C. SHELLARD 


Summary. A short description is given of Thom’s method of estimating the frequency of 
extreme winds. The method does not require data over a long period and is based on the 
expression G(v) for the mixed Fréchet type distribution of the extreme value v, where G(v) = 
be exp [—(o/B)—*) + pr exp [—(e/B)—* 5] in which the first term allows for extratropical 
storms and the second for tropical storms. pg and pv are the probabilities of an annual extreme 
being produced by an extratropical and a tropical storm respectively, and pg = 1 —py. The 
parameter § is related to the highest average monthly mean wind speed. 

An example is given of the detailed computation of the extreme-value wind distribution 
for Seawell, Barbados. For a selection of stations in the Caribbean, tables are given indicating 
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‘fastest mile’ speeds and some corresponding maximum 3-second gust speeds, for various return 
periods. The British Standard Code of Practice on Wind Loading has adopted the maximum 
3-second gust speed at 10 metres with a return period of 50 years as the basic wind speed for 
design purposes and corresponding basic design speeds are given for various Caribbean 
countries. 


Introduction. For the design of buildings and other structures the engineer 
requires information on the statistical probability of short-duration wind 
speeds. Ideally, rational design wind speeds should be based on the statistical 
analysis of suitable wind speed recordings covering periods of 20 years or 
more from a well-distributed network of meteorological stations. Unfortunately 
such data are available for only very few parts of the world. 


In the Commonwealth Caribbean, suitable wind records for the purpose 
are currently being obtained and, in some cases, have been available for 
quite a number of years. Until recently, however, few of these records have 
been subjected to routine measurement and tabulation and even if some of 
the longer records could be located and analysed it is fairly certain that no 
homogeneous set of records covering a sufficiently long period would result. 


Thus, if estimates of extreme wind speeds having stated probabilities are to 
be obtained for the Commonwealth Caribbean area now, some alternative 
procedure must be used. Fortunately such a procedure has been developed 
by Thom! and this paper is concerned with the application of his method 
to places in the Caribbean and the presentation of the results. 


The quasi-universal extreme wind distribution (Thom). Extreme- 
value wind distributions that have been established for areas with suitable 
long-period records, such as the U.S.A. (Thom?) and the U.K. (Shellard*:‘), 
are determined by two parameters: a shape parameter y and a scale para- 
meter 8. The shape parameter defines the shape of the frequency distribution 
and the scale parameter determines the wind speeds themselves. Thom has 
shown that for extratropical storms and thunderstorms in the U.S.A. the 
Fréchet distribution of the extreme value v, 


F(v) = exp[—(2/B)~], aren | ) 


has a shape parameter y which approaches a value of g as record length 
increases. He also found that the scale parameter 6 is well correlated with the 
highest average monthly mean wind speed. Hence a method of approximating 
the extreme wind distribution in extratropical storm areas was obtained. 
When the limited data available for tropical storm areas were examined 
Thom found that the scale parameter could be estimated in the same way 
as for extratropical storms but that the shape parameter was only about half 
that found for extratropical storms. This result was based on data from nine 
stations, only three of which, Miami (Florida), San Juan (Puerto Rico) and 
Hong Kong, had records covering more than 20 years. Nevertheless, Thom 
concluded, very reasonably, that since in the southern North Atlantic and the 
Caribbean some of the annual extreme wind speeds are due to extratropical 
storms (or thunderstorms) the annual extreme population will be a mixture 
of extremes from tropical and extratropical storms, giving a mixed extreme- 
value distribution function : 


G(v) = peFr(v) + prF r(v), are 
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where Fg(v) = exp[—(v/8)-*), ++ (3) 
Fr(v) = exp[—(v/8)-**] oo 
and pe = 1—p7, see (5) 


where pz and pr are the probability of an annual extreme wind being produced 
by an extratropical and a tropical storm, respectively. He went on to show 
that pr could be estimated from the mean number of tropical storm passages 
per year through the area concerned and presented a map of this parameter 
(Figure 3 of Thom’s paper'), covering the southern North Atlantic and the 
Caribbean and based on storm tracks for the years 1901-63. 


The relationships between 8 and vm, the highest average monthly mean 
wind speed, and between f7 and f, the mean number of tropical storms per 
year, established by Thom are : 


8 = (320°5 tm + 248°7)*—15°7, - ++ (6) 
br = 1/[1 + 99 exp (—3:0f )]. -++ (7) 


Application to stations in the Commonwealth Caribbean. Using 
Thom’s results, summarized above, it is clear that an extreme wind distribution 
can be estimated for any station for which vy, can be computed from available 
observations of wind speed, or reasonably estimated. Values are given in 
Table I for a number of stations in the Caribbean, together with the period 
of records used and the corresponding values of 8. Also given are the 
appropriate values of f taken from Thom’s map and the corresponding values 
of pr. San Juan was included in Table I so that extremes estimated by this 
method might be compared with those computed in the more :sual way 
(analysis of the annual extremes over a long period) and given by Thom in 
Table 2 of his 1968 paper.? 


TABLE I-—DATA REQUIRED FOR ESTIMATION OF EXTREME WIND DISTRIBUTION 
FOR CARIBBEAN STATIONS 
Station Um Period 8 
mile|h 
San Juan, Puerto Rico 12°9 1940-55 50°5 
Palisadoes, Jamaica 13°4 1950-62 51°7 
Coolidge, Antigua 15'0* 1941-48 55°4 
Sea 16°5 1954 58°7 


1g°0 1954-60 50°7 
lad I 1954-60 37°5 
East coast, Trinidad 12'0* 48°3 
Crown Point, Tobago 13'0* 50°7 


* Estimated. 


In three cases the values of vm, were estimated. The only average wind 
speed datum readily available for Coolidge, Antigua, was an annual average 
of 13:4 mile/h over the years 1941-48, and 15-0 mile/h is therefore a con- 
servative estimate of the highest average monthly mean wind speed. The 
mean speed given for Piarco was considered to be unrepresentative of Trinidad 
for the purpose of estimating a design wind speed since Piarco is rather 
sheltered and well inland in relation to the prevailing winds. A mean speed 
some 50 per cent higher was therefore estimated as representative of the east 
coast and hence of more exposed areas generally in the island, bearing in 
mind that winds associated with a tropical storm may blow from any direction. 
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The mean speed given for Crown Point was based on the analysis of 6 months 
of concurrent hourly wind speeds from Crown Point and Piarco, giving a 
mean ratio of 1-65 between the two stations. 

An example of the detailed computation of a mixed extreme-value wind 
distribution is given by Thom, but another is given in Table II using the 
data for Seawell, given in Table I. As Thom indicated, the distribution G(v) 
must be computed for a series of values of v and the speed that will be 
exceeded with any desired probability must then be obtained by interpolation. 
To compute F it is convenient to rearrange equation (1) in the form 


In In (1/F) = —y In (0/8), cee 
giving 1/F = exp exp [—y In (2/8)]. 


Also, since, for Seawell, 8 = 58-7 and pr = 0°13, we have to evaluate the 
mixed distribution 


G(v) =0'87 exp [—(0/58-7)-*] +0°13 exp [—(2/58-7)-**], +++ (10) 


and this is done in Table II for a series of values of v. 


TABLE II—COMPUTATION OF MIXED EXTREME WIND DISTRIBUTION FOR SEAWELL, 
BARBADOS 
Se In(v/58-7) —g In(v/58-7) —4:5In(o/58-7) 1/Fe 1/Fr O87FeE o13F7 G 
mile! 
3° —o-1602 1°4418 07209 «668-7 -7°815 Ss 013, 0°017 
00219 —0'1971 —0°0985 2273 2°475 0-383 0053 0°436 
0°81 


80 03097 — 27873 —1°9937 1-063 1-282 O'101 0-919 
100 0°5327 47943 —2°3971 1° 1095 0863 0118 0-981 
120 0°7149 — 6-434! —3°2171 1002 1°04! 0°869 0-125 0°994 
The values of v and G may be plotted on extreme-value probability paper 
as in Figure 1, from which the speed corresponding to any desired cumulative 
probability (or return period) may be read off. For example, the speed 


0999 ————-——-+ 
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having a cumulative probability of G = 0-98 is 100 mile/h. In other words 
a wind speed of 100 mile/h at Seawell is likely to be exceeded on the average 
only once in 50 years. 


Similar computations have been carried out for all the places listed in 
Table I with the results set out in Table III. It is important to emphasize 
here that these speeds are ‘fastest mile’ speeds, i.e. they are estimates of the 
fastest mile of wind. This is because Thom’s distribution is based on observa- 
tions made in the U.S.A., where the extreme wind speed variable that is 
most widely available is the fastest mile. The values given in brackets on the 
first line of Table III are the Puerto Rico fastest mile values published by 
Thom? and mentioned earlier. It will be noted that the agreement between 
the two sets of values is quite good. 


TABLE III—FASTEST MILE SPEEDS* FOR RETURN PERIODS OF 10, 20, 25, 50, 100 
AND 200 YEARS 


Return period (years) 
2 te) 100 


mile/h 
83 (80) 94(95) 105 (110) 


71 7 83 105 
78 gl 113 
79 gl 110 
67 77 94 
48 56 69 77 
63 72 86 93 
67 77 93 103 


* The speeds given in mile/h can be converted to knots by multiplying by 0-868, and to metres/ 
seconds by multiplying by 0-447. 
Figures in brackets are the Puerto Rico fastest mile values published by Thom.? 


In the West Indies the fastest mile of wind is not normally measured and 
the British practice of extracting the highest hourly mean speed and the 
highest gust speed for each day is generally followed. The anemographs in 
most general use have response characteristics such that the maximum gust 
speed is equivalent to a speed averaged over a period of approximately 3 
seconds. The fastest mile speeds given in Table III may be converted into 
their equivalent 3-second gust speeds, by making use of experimental measure- 
ments that are available on the relationship between wind speeds averaged 
over different periods of time. Such measurements, over open level country, 
have been made, for example, at Cardington, England (Durst®), and at Sale, 
Australia (Deacon®), the two sets of results being in excellent agreement. 
By making use of these and remembering that v mile/h is equivalent to a 
speed averaged over 60/v minutes (or 3600/v seconds), a simple relationship 
can be deduced between the fastest mile speed (v) and the corresponding 
3-second gust speed (g). This relationship is : 


g,= 1:09 v + 8 mile/h. ste (a8) 
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Hence the maximum (3-second) gust speeds corresponding to the fastest 
mile speeds of Table III may be calculated and these are set out in Table IV, 
which thus gives estimated maximum gust speeds for places in the Common- 
wealth Caribbean for return periods of 10, 20, 50 and 100 years. 


TABLE IV-—-MAXIMUM GUST SPEEDS FOR RETURN PERIODS OF I0, 25, 50 AND 
I00 YEARS 
Return period (years) 
20 50 100 
mile{h 
Palisadoes, Jamaica 94 110 123 
Coolidge, Antigua ¢ 104 120 132 
Seawell, Barbados 104 117 128 
Pearls, Grenada 89 101 111 
Piarco, Trinidad 67 76 83 
East coast, Trinidad 83 93 102 
Crown Point, Tobago 89 101 110 


It may be noted that the British Standard Code of Practice on Wind 
Loading (British Standard Code of Practice 3, Chapter V, Part II, 1970) has 
adopted the maximum 3-second gust speed at 10 metres (33 ft) above the 
ground, with return period of 50 years, as the basic wind speed for design 
purposes. Thus, corresponding basic design wind speeds for some Common- 
wealth Caribbean countries, based on the figures in column 4 of Table IV, 
might be suggested for adoption in any new code of practice for use in the 
area, as indicated in Table V. 


TABLE V—SUGGESTED BASIC DESIGN WIND SPEEDS FOR SOME COMMONWEALTH 
CARIBBEAN COUNTRIES 
mile/h 
Jamaica 
Leeward and British Virgin Islands 120 
Barbados, St Vincent, St Lucia 
Grenada, Tobago 100 
Trinidad go 
Guyana 50 
In the last line of this table a figure of 50 mile/h has been inserted for 
Guyana, which is not affected by tropical storms. This is based on the fact 
that in a record commencing in 1927, the highest gust speed recorded at 
Georgetown (Botanic Gardens) was one of 43 mile/h. 


Finally, in view of the fact that wind speeds of 150 to 160 mile/h have been 
measured in well-developed hurricanes, and these are means over periods 
of a minute or longer, it might be thought that the speeds in Table V are on 
the low side. However, it must be pointed out that the return period of a 
well-developed hurricane in any particular place is probably much greater 
than 50 years. Moreover, an engineer who wishes to design a structure to 
withstand the full force of a major hurricane may be expected to take this 
into account, and to use a design wind speed of about 160 mile/h or more, 
with a return period of perhaps 1000 years, depending on an assessment of 
the risk of failure in relation to the purpose, probable lifetime and cost of the 
structure concerned. 


Acknowledgement. This paper is published by kind permission of the 
Secretary-General of the World Meteorological Organization, having been 
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551.509.325 551.591 
VISIBILITY DETERIORATIONS DURING WINTER MORNINGS 


By W. E. SAUNDERS 

Summary. A study of visibility deteriorations on clear winter mornings at a number of 
aerodromes in eastern England shows that the probability of deterioration varies markedly 
with the position of a site in relation to smoke sources and low-lying areas. Wind direction 
is therefore of great importance. The likelihood of deterioration is shown to decrease with 
increasing wind speed. It is also shown that on the majority of occasions when the lowest 
visibility reached during the morning did not exceed 2 km the wind was calm or the direction 
was from smoky or damp, low-lying areas. The time of reaching the lowest visibility is shown 
to vary from station to station. 


It is well known to forecasters that following a clear night, which has been 
free from fog, visibility frequently decreases during the morning, especially 
during the winter. This may be due to industrial or domestic smoke, or in 
some locations to the drift of mist or fog from low-lying or fen areas. In either 
case, the diurnal increase in wind speed and turbulence after sunrise con- 
tributes to the effect. 


The records for 10 aerodromes between Northumberland and the Cambridge 
area (see Figure 1) for the months November to February for five winters 
(November 1965 to February 1970), were examined in order to assess the 
incidence of these morning deteriorations. Cases included were those in 
which at 04 GMT there was no fog and the cloud amount did not exceed 3/8, 
excluding broken cirrus. The initial time of 04 Gmt was selected because this 
is the last observation available before forecasts require to be issued for many 
day-time activities. The sunrise times vary between about 0645 and 0830 GMT. 


For the purpose of this study, a deterioration was regarded as having 
occurred if the visibility decreased from within one of the following ranges 
at 04 GMT into any lower range within the period 06-13 GMT: 

1:0 — 1:9 km 
2:0 —- 3:9 km 
4:0 —- 8-0 km 
over 8 km 

Cases where a visibility deterioration was due to precipitation were 
omitted. 
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The overall position regarding deteriorations is given in Table I. This 
shows considerable consistency north of about York, with deteriorations 
on around 30 per cent of the clear mornings. The probability of deterioration 
increases in south Yorkshire, because of the proximity of Church Fenton and 
Finningley to the smoke sources of the West Riding and urban areas further 
south, and of the low-lying ground near these aerodromes. However, the 
difference between the frequency of deteriorations at stations Acklington 
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TABLE I—INCIDENCE OF DETERIORATIONS IN VISIBILITY BETWEEN 06 AND 13 GMT 
ON MORNINGS WHICH WERE CLEAR*® AT 04 GMT, NOVEMBER 1965 TO FEBRUARY 
1970 

Station Number of days Incidence of deteriorations 
clear at 04 GMT Number of days Percentage 

Acklington 164 44 27 
Leeming 173 51 29 
Topcliftet 133 38 29 
Linton-on-Ouset 105 ‘ 31 
Church Fentontt 96 8 
Finningley 18g be 
Manby 169 c 53 
Strubbyt 109 38 
Cranwellt 109 37 
Oakingtont 131 40 

* No fog and cloud not exceeding 3/8, excluding broken cirrus. 

+ These stations normally report on Mondays to Fridays only. 

¢t At Church Fenton 06 oat was taken as the initial time instead of 04 Gar. 


to Linton-on-Ouse and the frequency at Church Fenton and Finningley is 
also partly due to the original selection of clear nights. Cloudless weather 
during the night and early morning is most likely at Acklington and in north 
Yorkshire with gradient winds from between west and north, and the figures 
reflect lack of smoke sources within this quadrant for stations down to Linton. 
Also, of course, these airstreams have crossed the Cheviot and Pennine hills, 
and are relatively dry. The stations south of Finningley are affected by a 
variety of smoke sources, and by drift from low-lying ground, but the likelihood 
of deterioration is rather less than in south Yorkshire. 

Table II shows the effect of wind speed on the incidence of deteriorations. 
In general the probability of deterioration decreases with increasing wind 
speed. 

The effect of wind direction is shown in Table III. For each station, 
omitting the days of calm, the occasions were separated according to mean 
wind direction in the period from 06 Gar to the time of lowest visibility; those 
with wind directions from smoky or low-lying areas were separated from 
those with ‘clean’ wind directions. An attempt to separate the smoky occasions 
from those with drift from low-lying areas by use of a relative-humidity 
criterion was not satisfactory. Table III shows that wind direction is of great 
importance at Vale of York stations, especially at Leeming. The. results 
imply that the smoke trails which affect that aerodrome are very sharply 
defined. South of York the differences due to wind direction become less 
pronounced. 


In Table IV the lowest visibility in the period 06-13 cmt has been tabulated 
against wind speed, with the mean directions grouped as in Table III. This 
emphasizes the importance of smoke pollution, or in some cases drift from 
low-lying ground, in the formation of fog on winter mornings. Out of 129 
occasions of morning fog formation, the mean wind direction was from smoky 
or low-lying directions, or calm, on 115 occasions, and from ‘clean’ directions 
on 14 occasions. Similarly, out of 116 occasions when the lowest visibility 
was 1:0-1'9 km, the mean wind direction was from smoky or low-lying areas, 
or calm, on 100 occasions, and from ‘clean’ directions on 16. Even allowing 
for the fact that there were more total occasions with winds from smoky or 
low-lying areas, or calm (857), than from ‘clean’ areas (521), there is the 
strong implication that pollution is a major factor in deteriorations. 
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There is no doubt that the wind directions which have been taken as 
basically ‘clean’ in this investigation often carry some pollution from minor 
sources, and it is important in forecasting to identify and allow for these 
sources. It is also true that topographical effects sometimes lead to smoke 
arriving at an aerodrome from a direction which is normally ‘clean’. For 
example, at Topcliffe the wind direction for arrival of West Riding smoke 
is normally 180—220°, but with light winds under a winter inversion the smoke 
may reach the station from the south-east, having been confined to the east 
side of the Vale of York by high ground further east. 


In Table V, the time of reaching the lowest visibility within the period 
06-13 GmT has been related to certain wind directions for a selection of 
stations. For Leeming, with wind directions which bring West Riding smoke, 
there are maxima at 07-08 and 11-12 cmt. For Finningley, with winds from 
smoky or low-lying directions, there is a pronounced maximum of lowest 
visibility between 08 and 10 cmt. For Manby, the times have been given 
separately for the direction bands 010-170° (wind off the sea, or from damp, 


TABLE V—VARIATION OF TIME OF INCIDENCE OF LOWEST VISIBILITY IN PERIOD 
06 To 13 GMT WITH CERTAIN WIND DIRECTIONS ON MORNINGS WHICH WERE 
CLEAR AT 04 GMT, NOVEMBER 1965 TO JANUARY 1970 

Time of lowest visibility, Gut 
Station Mean wind 06-07 07-08 08-09 09-10 10-11 II-I2 12-13 
direction 
degrees true Incidence of lowest visibility (number of days) 
Leeming 140-190 2 6 3 2 I 8 


5 

9 

low-lying ground), and from 280-360° (predominantly smoky). This shows 
that with direction 010-170° the lowest visibility is reached in the early 
morning, but that with 280-360° there are maxima at later times, including 
one as late as 12-13 GMT. This emphasizes that deteriorations due to smoke 
from remote sources tend to occur in the late morning. 


Finningley 180-330 22 19 4! 47 5 I 
I — 
2 


Manby { 010-170 4 2 4 


280-360 7 11 12 


Acknowledgement. The author is indebted to outstation staff for their 
considered opinions on the smoke sources, and for extracting the data. 


551-554°551-550:551-453 
SAND STREETS 


By J. D. HASTINGS 


Summary. During a flight over a ‘ghibli’ between El Adem (Libya) and the coast on the 
morning of 14 February 1969, it was noticed that lines were visible below in the rising sand 
in a number of places. The lines were analogous to cloud streets aligned in the direction of 
the winds. 

It is suggested that what had been observed was the process of longitudinal dune formation 
in action. Study of the upper air data for that day shows that the diameter of helical eddies 
would have been about 300 m and the spacing between them 600 m to goo m, that is to say 
about half a nautical mile. 


A Hercules aircraft of Air Support Command took off from RAF El Adem 
at about 0930 GmT on 14 February 1969 in conditions known locally as 
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‘ghibli’, ie. the south-south-easterly winds were gusting at times to 40 kt* 
and visibility was reduced to a few hundred metres by rising sand. Not long 
after take-off when permission was given to release seat belts and when as 
far as could be estimated the aircraft was approaching the Jebel Akhdar 
to the south of Derna at a height of 10 000 ft, it was noticed on looking out 
of the port-side window that lines showed up quite clearly in the rising sand 
below; they occurred in a number of places and appeared to be aligned 
along the direction of the winds. Unfortunately, external vision is very 
limited in a Hercules aircraft except from the flight deck so it was not possible 
to form any idea of how widespread the phenomenon was, nor was there a 
camera readily available to photograph it. 


Initial reaction was to speculate whether the phenomenon showed up on 
pictures by automatic picture transmission (APT) but although an examination 
was made at E] Adem of such pictures of the area on occasions of ghibli going 
back almost a year, it proved quite fruitless. At the time it was not appreciated 
that the distance between the streets (approximately 0-5 to 1 km) would 
be so much less than the present resolution (4 km) of the satellite cameras. 
Later in the year while attempting to explain the phenomenon, not entirely 
satisfactorily, on the basis of ‘thermal streeting’ (Konrad!) the author’s 
attention was drawn to a paper by Hanna.? In it he states that satellites 
orbiting 200 to 300 km above the earth have confirmed the presence of 
longitudinal dune systems in nearly every large desert area of the world. 


Satellite pictures show the dunes aligned north-north-west to south-south- 
east over Libya and Egypt, turning north-east to south-west over Chad and 
the Sudan. This direction over Libya is in accord with statistics which show 
that for El] Adem 50 per cent of winds at 06 and 12 ont of force four or more 
are from 330/360° and 11 per cent from 150/180°. The spacing over the 
Fezzan is estimated to be 2 to 3 km but in general over the world the average 
is 1 to 2 km. 


Longitudinal vortices have been observed in the motion of neutral-lift 
balloons. Angell et alii* observed helical circulations of 600-m diameter over 
the Los Angeles basin and in 1968‘ noted counter-rotating helical circulations 
over the Idaho desert having a diameter of 1-5 to 2-5 km. 


Cloud-street spacings appear to range from o-5 km to as much as 30 km. 
Kuettner® made a thorough investigation into cloud streeting and found 
the average spacing to be between 5 and 10 km, a value usually twice to three 
times the depth of the convective layer. He gives the following conditions 
as conducive to streeting : 

(i) flat underlying terrain, 

(ii) little variation of wind direction with height, 

(iii) wind speed higher than normal, 

(iv) strong curvature of wind-speed profile, and 


(v) unstable lapse rate close to the surface (and an inversion at the top 
of the convective layer). 





* Conversion factors to metric units are: 1 knoto-5 m/s; 1 foot = 0-3048 m. 
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When these conditions are met the dominant form of motion in the 
convective layer is a longitudinal roll vortex approximately parallel to the 
mean flow and having a wavelength about twice the thickness of the layer. 


After considering alternative approaches to the problem of helical eddy 
formation Hanna believes Kuettner’s theory to be the most applicable to 
the boundary layer of desert areas. The process envisaged is shown in 
Figure 1. 


Cloud Cloud 


ia. 


D oe” 
tj WH Y=Z@#!?!H- 


FIGURE I-——-MECHANISM OF DUNE AND CLOUD FORMATION 


It was therefore decided to see how well Kuettner’s five conditions were 
satisfied on the occasion when the sand streets were noted. The wind profiles 
for El Adem for 05 Gmr and Tobruk for oo, 06 and 12 Gat are shown in 
Figure 2. For comparison the mean wind profile found by Kuettner for 
cloud streeting is included. The ascents from Tobruk for oo and 12 GMT 
are given in Figure 3. 
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FIGURE 2—-WIND PROFILES FOR TOBRUK AND EL ADEM, I4 FEBRUARY 1969 


Hanna has noted that the wind must pass over a minimum distance of 
flat surface before longitudinal vortices will develop and, using the time 
period T of the circulations within observed vortices of about 2000 s and a 
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FIGURE 3—ASCENT TEMPERATURES FROM TOBRUK FOR 00 AND I2 GMT, 
14 FEBRUARY 1969 


typical wind speed U of about 10 m/s, arrives at a value for this minimum 
of J x U = 20 km approximately. Although there is some uncertainty as 
to the exact location where the sand streets were seen it is thought that 
condition (i) as well as the requirement of a ‘fetch’ of more than 20 km were 
both satisfied. 


From the wind-speed profiles it will be seen that at both El Adem and 
Tobruk at 06 cnr in particular, there is only a small variation of direction 
with height in the lowest kilometre. By 12 Gar although direction is still 
unchanged in the lowest half kilometre a veer of 30° has taken place at 1 km 
and above. Condition (ii) is however met within the convective layer. 


Surface winds at El Adem and Tobruk were 20 kt at 05/06 cmt and had 
increased to 25/30 kt with gusts to 40 kt from 160° by 09 cmt. Condition 
(iii) is evidently satisfied. 


It will also be seen from the wind-speed profiles that a maximum speed 
was reached at a height of about 300/450 m followed by a decrease of speed 
to about 1500 m. This effect is still noticeable on the 12 car profile, although 
much less pronounced, and the decrease above 1 km no longer obtains. 
Condition (iv) was certainly met at 05/06 cmt and most probably at the 
time of sighting the streets, but only to a lesser extent by 12 Gr. 


Finally, the lapse rate near the surface, condition (v) : by interpolation 
between the two ascents it is evident that, taking the og GT surface temperature 
of 20°C, a marked inversion from about 300 m to 1 km must have existed 
and that therefore convection must have been restricted to the layers below it. 


Applying the empirical relationship of the spacing being two to three 
times the depth of the convective layer gives spacing values of 600 to goo m 
or about half a nautical mile. These figures correspond to the lowest of the 
range of spacings for cloud streets or longitudinal dunes. 


These distances are well below the resolution possible at present on APT 
pictures but not on high-level aerial photographs such as made from Gemini 
spacecraft. The National Aeronautical Space Administration (NASA) Manned 
Spacecraft Centre at Houston, Texas, provided three colour transparencies 
of the south-west of Libya which were the only photographic coverage 
available of the country. An approach was then made to the Joint Air 
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Reconnaissance Interpretation Centre (JARIC) but again, unfortunately, 
without success. Finally however, a series of photographs covering the area 
was obtained from the Survey Liaison Staff (UK) Royal Engineers (SLS 
(UK)RE) from which it was evident both to skilled and unskilled eyes that 
no permanent dune formations could be seen at least at the time that the 
photographs were taken. 


From the most recent maps and from personal knowledge of the area it is 
apparent that while there is sufficient sand to produce serious deteriorations 
of visibility with strong winds it is not enough for permanent dunes to form. 
Any which did during the period of a ghibli could be expected to be of little 
vertical extent and of a temporary nature, being easily destroyed subsequently 
by winds blowing from another direction. Nevertheless, it has been shown 
that the conditions for longitudinal dune formation as postulated by Hanna 
were satisfied on the morning of 14 February 1969 and it is suggested that what 
appeared to be sand streets was, in all probability, this process as it would 
appear from above. 


The author would like to acknowledge with thanks the transparencies 
provided by NASA and the help given by both JARIC and SLS(UK)RE. 
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REVIEW 


Atmospheric optics (translated from the Russian), edited by Nikolai B. Divari. 
270 mm X 210 mm, pp. 191, ilus., Consultants Bureau, London, Donnington 
House, 30 Norfolk Street, London WCa2, 1970. Price: $27-50. 


This book contains the proceedings of conferences on atmospheric optics 
sponsored by the committee on the optical instabliity of the atmosphere of 
the Astronomical Council of the Academy of Sciences of the U.S.S.R. The 
meetings were held in Pulkovo in November and December 1965 and in 
November 1966. 


There are 31 short papers, about 10 of which are of mainly astronomical 
interest. The remainder deal with atmospheric transmission coefficients, 
atmospheric scattering, spectral brightness of the sky, polarization and 
twilight phenomena, including the evaluation of dust in the upper atmosphere 
by the twilight method. 

J- PATON 
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OBITUARIES 
It is with regret that we record the death on 10 February 1971, of Mr H. J. 
Groom, Experimental Officer, Met O 8, and the death on 5 February 1971, 
of Mr J. W. Armstrong, Senior Scientific Assistant, Acklington. 


OFFICIAL PUBLICATIONS 
Geophysical Memoirs 
No. 113. An observational study of the meridional flux of energy and angular 


momentum in the troposphere and lower stratosphere at latitude 30°N using 
1958 IGY data. By A. E. Parker. 


This memoir discusses the results of the calculations made from information 
amassed during the International Geophysical Year, 1958, concerning the 
transport of momentum, heat and energy northwards across latitude 30°N. 
Four months — March, June, September and December — were chosen as 
being the most representative for the year. 


Scientific Paper 
No. 31. The three-dimensional analysis of meteorological data. By R. Dixon, 
B.Sc. and E. A. Spackman, M.Sc. 


Before meteorological observations can be used by a computer for the 
purpose of forecasting the weather, the raw observations have to be processed 
within the computer so that they can serve as initial data for the mathematical 
model of the atmosphere which actually does the forecasting. The advent 
of modern observational devices such as satellites and long-life free drifting 
balloons has greatly complicated the process of data analysis by computer. 
This paper presents a possible method for effecting this data analysis by 
fitting a high-power polynomial to all the observations from within a large 
volume of the atmosphere. 


The following publication has recently been issued : Monthly charts of dew- 
point temperatures over the Indian Ocean. 


This publication consists of charts of the mean and the 95 and 5 percentile 
surface dew-point temperatures over the Indian Ocean for each month of 
the year. It is an addition to Monthly meteorological charts of the Indian Ocean 
(HMSO, 1949). These dew-point temperature charts were prepared from 
punched cards of ships’ observations, mainly British, for the period April 
1953 to December 1961. This representation of humidity will be useful to 
shippers and packers who are concerned with problems of cargo ventilation 
in the Indian Ocean area. 
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‘RADIO SONDE 


* Meteorological Transmitter 


The WB Radio Sonde is essential for high altitude weather recording 

(up to 66,000ft.), and is available with parachute, radar reflector and 

battery, or as a single unit, complete with met. elements. 

For full specification of the WB Radio Sonde—which is used by the U.K. 
Meteorological Office, and many overseas Governments —please write or telephone 








MANSFIELD 
NOTTS 


FLECTRIGAL RADIO CO. LTD. nnn. 
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